Doped Kagome System as Exotic Superconductor 
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A Chern-Simons theory for the doped spin-1/2 kagome system is constructed, from which it is 
shown that the system is an exotic superconductor that breaks time-reversal symmetry. It is also 
shown that the system carries minimal vortices of flux hc/Ae (as opposed to the usual hc/2e in con- 
ventional superconductors) and contains fractional quasiparticles (including fermionic quasiparticles 
with semionic mutual statistics and spin-1/2 quasiparticles with bosonic self-statistics) in addition 
to the usual spin-1/2 fermionic Bougoliubov quasiparticle. Two Chern-Simons theories — one with 
an auxiliary gauge field kept and one with the auxiliary field and a redundant matter field directly 
eliminated — are presented and shown to be consistent with each other. 



I. INTRODUCTION 

The "perfect" spin-1/2 kagome lattice, realized re- 
cently in Herbertsmithite ZnCu3(OH)6Cl2r'i^i^ has pro- 
duced great enthusiasm in both the experimental and the 
theoretical condensed matter community. Experimen- 
tally, the antiferromagnetic exchange is found to be J « 
190 K, and yet no magnetic ordering is observed down 
to a temperature of 50 mK.— Theoretically, with nearest- 
neighbor Heisenberg antiferromagnetic interaction, sev- 
eral possible ground states have been proposed, includ- 
ing the valance bond solid (VBS) states^ and the Dirac 
spin liquid (DSL) statej^ii while results from exact diag- 
onalization (ED)^ remains inconclusive as to which state 
is preferred. 

So far both the experimental and theoretical studies 
have been focused on the half- filling (i.e., undopcd) case. 
In this paper, we investigate the situation in which the 
kagome system is doped, which could in principle be re- 
alized by substituting CI with S. We shall take the DSL 
state, which at low energy is described by spin-1/2 Dirac 
fermions (spinons) coupled to an emergent internal gauge 
field, as our starting point. Naively, one might expect 
the system to be a Fermi liquid with small Fermi pock- 
ets opening up at the spinon Dirac nodes. However, 
since the system contains an emergent internal gauge field 
a^, filled Landau Levels (LLs) can spontaneously form. 
When the flux quanta of this emergent field is equal to 
half of the doping density, the resulting LL state is en- 
ergetically favorable, (the formation of filled LLs, as in- 
duced by the internal gauge flux, has also been proposed 
in the case when an external magnetic fleld is applied 
to the undoped spin-1/2 kagome system) i^ Furthermore, 
the strength of this a field and the doping density can 
co-fluctuate smoothly across space, resulting in a gap- 
less excitation in density mode. Since the gapless den- 
sity mode is the only gapless excitation, the LL state is 
actually a superconducting state. This provides an un- 
conventional superconducting mechanism which results 
in a time-reversal symmetry breaking superconductor. 

As typical for a superconductor, the state we proposed 
also supports electromagnetically (EM) charged vortices. 
In additional, since there are multiple species of emergent 



spinons and holons, the system also contains EM-ncutral 
topological excitations that arc analogous to quasipar- 
ticles in quantum Hall systems. To describe the super- 
conducting state, the EM-charged vortices, and the EM- 
neutral quasiparticles in a unified framework, we start 
with the t-J model and the DSL ansatz and construct 
a Chern-Simons theory, well-known from the study of 
quantum Hall systems, for this system. 



In our scenario, the low-energy effective theory con- 
tains four species of emergent holons, each carry a charge 
e. All four species are tied together by the emergent 
gauge field a^ . Consequently, the flux of a minimal vor- 
tex in this superconductor is found to be hc/Ae, as op- 
posed to the usual hc/2e in conventional superconduc- 
tor. Furthermore, the quasiparticles in this scenario are 
shown to exhibit fractional statistics. In particular, there 
are fermionic quasiparticles with semionic mutual statis- 
tics and bosonic quasiparticles carrying spin 1/2. 



This paper is organized as follows: In Sect. [TTl we 
derive the Chern-Simons theory starting with the t-J 
model and motivate the necessity of such an "uncon- 
ventional" formation for superconductivity. In Sect. IIIIl 
the existence of superconductivity is first explained intu- 
itively, and then confirmed by a more rigorous derivation. 
The physical vortices are then discussed, with the hc/Ae 
magnetic flux explained both intuitively and mathemat- 
ically. In Sect. IIVI the EM-neutral quasiparticles are in- 
troduced and their statistics are derived. The discussion 
on these quasiparticles continue into Sect. |V] in which 
their quantum numbers are analysed. In Sect. IVIl an 
alternative formulation of the Chern-Simons theory is 
presented, in which the auxiliary gauge field a^ and a re- 
dundant matter field are eliminated directly, and the re- 
sults obtained are shown to be consistent with that of the 
previous sections. The paper concludes with Sect. IVIII 



II. FROM t-J HAMILTONIAN TO 
CHERN-SIMONS THEORY 

The starting point of our model for the doped kagome 
system is the t-J Hamiltonian: 

HtJ^^Jh>^■ Sj - -n.rij j -t{c\^Cj„ + h.c/j , (1) 

where cj^. and Cjo- are projected electron operators that 
forbid double occupation, and that J > 0. Throughout 
this paper we shall assume that t > and that the system 
is hole doped. For t < 0, our results can be translated to 
an electron-doped system upon applying a particle-hole 
transformation. 

Using the C/(l) slave-boson formulation^i^ we intro- 
duce spinon (fermion of charge and spin 1/2, represent- 
ing singly occupied sites) operators ji„ and holon (boson 
of charge -\-e and spin 0, representing empty sites) oper- 
ators hi such that c^^ ~ fJa^i' ^^^ ^-PPly the Hubbard- 
Stratonovich transformation. This yields the following 
partition function: 



I DfDpDhDh*D\DxDAcxpi- I drLi 
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in which the mean-field conditions are given by Xij = 

Y.a(fLfj<y) and A^j = {f^-ffji - /a/iT>; 

Assuming mean-field ansatzes in which Aij = and 
Xij = X6~*"'^i and rewriting A^ = ckq, we arrive at the 
following mean-field Hamiltonian: 

+ E '*l(*"o - ^^B)hi - tx^ie'""' hlhj + h.c.) . 

(4) 
Observe that an internal gauge field a'^ emerges nat- 
urally from this formulation. Its space components aij 
arise from the phases of Xij i while its time component oq 
arises from enforcing the occupation constraint: 
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FIG. 1: (a) The kagome lattice with the DSL ansatz. The 
dashed lines correspond to bonds with t = — 1 while unbroken 
lines correspond to bonds with t = 1. ri and r2 are the 
primitive vectors of the doubled unit cell, (b) The original 
Brillouin zone (bounded by unbroken lines) and the reduced 
Brillouin zone (bounded by broken lines) of the DSL ansatz. 
The dots indicate locations of the Dirac nodes at half-filling 
while the crosses indicate locations of the quadratic minima of 
the lowest band, ki and k2 are the reciprocal lattice vectors 
of the reduced Brillouin zone. 



From Eq. 31 it can be seen that the holons and spinous 
arc not directly coupled with each other at the mean- 
field level — they arc correlated only through the common 
gauge field a^. Consequently, if we treat a^ at the mean- 
field level, the spinon spectra and the holon spectra will 
decouple, and up to an overall energy scale both will be 
described by the same tight-binding Hamiltonian. 

By gauge invariance, a mean-field ansatz for a'^ is 
uniquely specified by the amount of fluxes through 
the triangles and the hexagons of the kagome lat- 
tice. In particular, the DSL state is characterized by 
zero flux through the triangles and tt flux through the 
hexagonsi^ii By picking an appropriate gauge, the DSL 
state can be described by a tight-binding Hamiltonian 
with doubled unit cell, in which each nearest-neighbor 
hopping is real, has the same magnitude, but varies in 
For the precise pattern see Fig. 



sign, ]^or tne precise pattern see i*ig. |l(a) This tight- 
binding Hamiltonian produces six bands, whose disper- 
sions are, in units where the magnitude of the hopping 
parameter is set to 1, 



^-top — 2 



(doubly degenerate) 



(6) 
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1 ± \/ 3 =F v2\/3 — cos 2kj: + 2 cos k^ cos v3fcy 

(7) 



At any k-point, E^,+ < S_,_ < S+,_ < E+^+ < Stop- 
These tight-binding bands have the following features 
that will be important for our purposes: (1) four degen- 
erate shallow quadratic band bottoms in the flrst (low- 
est) band i?_,+ ; and (2) two degenerate Dirac nodes 
where the third band (£'+__) and the fourth band (i?+,+) 
touches. See Fig.|l(b) and Fig. [5] for illustrations. 



Kh^+fy^1+fya^^■ 



(5) 



Now suppose the doped kagome system is described 
by the DSL ansatz as in the undopcd case, and that the 
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FIG. 2: The band structure of the kagome lattice with the 
DSL ansatz (a) plotted along the line kx = and (b) of the 
lowest band plotted along the line ky = 0. Note that the top 
band in (a) is twofold degenerate. 



doping is x per site. Then each doubled unit cell will 
contain 6x holons and 3 — 3x spinous per spin. By Fermi 
statistics, the spinous will fill the lowest 3 — 3a; bands 
and thus can be described by auti-spinon pockets at each 
Dirac node. Similarly, by Bose statistics the holons will 
condense at each quadratic band bottom. This state shall 
be referred to as the Fermi-pocket (FP) state. 

However, the FP state is not the only possibility. In 
particular, an additional amount of uniform a field can 
be spontaneously generated to produce LLs in both the 
holon and spinon sector. The resulting state shall be re- 
ferred to as the LL state. In the absence of holons (i.e., 
at half filling), both mean- field calculation and projection 
wavefunction study indicate that the LL state is energet- 
ically favored over the FP statci^ Since the spinon bands 
are linear near half-filling while the lowest holon band is 
quadratic near its bottoms, at the mean-field level the 
energy gain from the spinon sector (which scales as 3/2 
power of the a field strength) will be larger than the 
energy cost in the holon sector (which scales as square 
of the a field strength) at low doping. Therefore, even 
after the holons are taken into account, the LL state is 
expected to have a lower energy than the FP state. 

Furthermore, from mean-field it can be seen that the 
energy gain will be maximal when the a field is adjusted 
such that the zeroth spinon LLs are exactly empty. Since 
each flux quanta of the a field corresponds to one state in 
each LL, and that each anti-spinon pocket contains 3a;/2 
states for a doping of x per site, the fiux must be 3a; flux 
quanta per doubled unit cell for the zeroth spinon LLs to 
be empty. 

As for the holon sector, there are 6x holons per dou- 
bled unit cell or equivalently 3a;/2 holons per band bot- 
tom. Since the holon carries the electric charge and are 
hence are mutually repulsive, one may expect them to flU 
the four band bottoms synmretrically. In such case the 
flrst LL of each of the holon band bottom would be ex- 
actly half-filled, which implies that the holons would form 
four Laughlin ly — 1/2 quantum Hall states. Since the 
Laughlin z/ = 1/2 state is gapped and incompressible, this 
symmetric scenario should be energetically favorableii^ 

From the physical arguments given above, it can be 



seen that the effective description of this system is anal- 
ogous to that of a (mulit-layered) quantum Hall system, 
and thus may contain non-trivial topological orders, man- 
ifesting in, e.g., fractional quasiparticles with non-trivial 
statistics. In order to describe such system, we adopt a 
hydrodynamic approach well-known in the quantum Hall 
literature i^i^ii^ In this approach, a duality transforma- 
tion is applied, in which a gauge field is introduced to 
describe the current associated with a matter field, and 
which the two are related by: 



Zn 



(8) 



where J^ is the current of the matter field and a\ is the 
associated gauge field. Here fi, v, and A are spacetime 
indices that run from to 2, and e^^^ is the totally an- 
tisymmetric Levi-Civita symbol. 

In this formalism, a single-layer quantum Hall system 
of flUing fraction (a.k.a. Hall number) v = 1/m is de- 
scribed by the following effective Lagrangian: 



^ = - — ^''''^aMa.aA-— e'^^'^aM^-^A+^aMJ^ 



47r 



27r 



(9) 



where A^ is the external electromagnetic field and jy is 
the current density associated with particle-like excita- 
tions. The ". . ." represents terms with higher derivatives, 
and hence unimportant at low energies. In particular, at 
the lowest order in derivatives among the terms dropped 
is the "Maxwell term" : 
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1 
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{d^,a,-d,a;){d^a'' -d^a^). (10) 



The effective Lagrangian Eq. [9] can be understood by 
considering the equation of motion (EOM) with respect 
to the dual gauge field a^ . With a stationary quasiparti- 
cle at xq such that jy = ((5(a; — a;o),0, 0), the EOM reads, 
in the time-component: 



27r 



(11) 



which confirms that i/ indeed equals to the filling fraction 
27rJ°/(— ei?), and that jy = {6{x — Xo),0,0) is a source 
term for a quasiparticle having charge iv. In particular, 
a physical electron at Xq can be associated with jy = 
{S{x-xo),0,0) and £ = 1^-1. 

Since jy is a source of "charge" in a^, from the duality 
transformation Eq. [51 it can alternatively be viewed as a 
source of vortex in the matter field current J'' . 

The statistics of the quasiparticles can be deduced by 
integrating out the dual gauge field a^ in Eq. [9l from 
which we obtained the well-known Hopf term: 



C = nj^'v 



^fivxd" ^ ~:^ 
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(12) 



where j'' = 

that couple linearly to a^ 



e/2T:)e'"^^d''A + £j!^ is the sum of terms 



The statistical phase 9 when one quasiparticlc de- 
scribed by £ = -^1 winds around another with i = £2 
can then be computed by evaluating the quantum phase 
giS ^ e'-/'^', with ji" = tij!^^ + I2JV2 being the total 
current produced by both quasiparticles. This yieldsi^ 

In particular, for the statistical phase accumulated 
when an electron winds around a quasiparticlc of charge 
Iv to be a multiple of 27r, I must be an integer. This 
provides a quantization condition for the possible values 
ofl 

For an iV-layer quantum Hall system, Eq.[H] generalizes 
to: 



C 



+ (-laifj.Jv + • • ■ 
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+ {I ■ a,,)j(; + . . . 



(13) 



here a^ is the dual gauge field corresponding to the 






)' and 



matter field in the /-th layer, a^ = (a^, 
q = (gi, . . . , qn)'^ are A^-by-l vectors, £ = i^^i, . . . ,tN) 
is an iV-by-1 integer vector, and K = [Kjj] is an iV-by- 
N real symmetric matrix. On the second line of Eq. [T2] 
and henceforth, we adopt a condensed notation in which 
the boldface and dot-product always refer to the vector 
structure in the "layer" indices and never in the space- 
time indices. 

In the multi-layer case, assuming that dot K ^ 0, the 
procedure for integrating out the dual gauge fields can 
similarly be carried out, which yields: 



£' = TTifrh- 



kO- 



Q2 



(14) 



where j'' = -q(e/27r)e^''^9'M + ^j(t. The statistical 



phase 9 when one quasiparticlc described by .^ = ^1 winds 
around another with (. = (.2 can then be computed in 
a similar way as in the single-layer case, which yields 
9 = 2tt ij K^^i2- The information of quasiparticlc statis- 
tics is thus contained entirely in K^^. 

Except for the complication that there is both an ex- 
ternal EM field A^^ and an internal constraint gauge field 
a^ , the doped kagome system we proposed is completely 
analogous to a multi-layer quantum Hall system. We 
shall therefore construct a Chern-Simons theory similar 
to that of Eq. [14] by assigning a dual gauge field to each 
species of matter field. 

For the holon sector, wc can represent the holons at 
each of the four band bottoms by a dual gauge field bj 
(J = 1,2,3,4). Since the holons at each band bottom 
form a Laughlin v ~ 1/2 state, the total Hall number for 
the holon sector is ^ j vj = 2. For the spinon sector the 
situation is more subtle. Since the zeroth LL is empty 
and all the LLs below it arc fully filled at each Dirac node, 
wc may represent the spinous near each of the four Dirac 
nodes by a dual gauge field a^ (/ = 1,2,3,4) having Hall 
number v ~ —1. However, since a is internal the com- 
bined system of holons and spinous must be a neutral, 
which requires X^aii species ^^ = and hence in the spinon 
sector ^j lyj = —2. To circumvent this problem, we in- 
troduce two additional dual gauge fields ag and ag , each 
having Hall number v — +1. The two fields Og and Og 
can be thought of as arising from the physics of spinous 
near the band bottoms of the two spin species. In this 
setting, a^, . . . , 04 are expected to carry good spin and 
k quantum numbers fi^ while a^ and Og are expected to 
carry good spin quantum number only. Note also that 
a^, . . . , 04 possess an emergent SU{A) symmetry of spin 
and pseudo-spin (i.e., k-points). 

Assembling the different species, the low-energy effec- 
tive theory for the doped kagome system is given by the 
following Chern-Simons theory: 



1 '^ 1 '^ 9 1 / \ 

"" 7=1 ^ 7=5 "" J "" \ 7 J / 



|- ^ e^''^6,^9,AA + ( E ^^«^^ + E ^^^•^'^ ) 
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(15) 
(16) 



As before, the ". . ." denotes terms higher in derivatives, including first and foremost the Maxwell term analogous to 
Eq. [To] In the second line, we have combined the eleven gauge fields internal to the system into a column vector 



c'^ = (a'^;a^, 



J "-6 ' "1 ' 



b^. 



. ,64)'^- Note that unlike Eq. [T31 we have included the internal gauge field a'^ in c'*. This 



is because a^ is internal and can be spontaneously generated while the EM field in the usual quantum Hall case 
is external and fixed. This distinction is crucial, as will be evident soon. The "charge vector" q in this case is 



q = (0; 0, 0, 0, 0, 0, 0; 1, 1, 1, 1)^, and the K-matrix K takes the block form: 
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(17) 
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The three terms in Eq. [16] can be understood as fol- 
lows: the first term describes smooth internal dynamics 
of the system; the second term describes its response un- 
der an external EM field; and the third term describes 
the topological excitations of the system, which can be 
thought of as combinations of vortices in various matter- 
field components. As in Eq. [131 ^ must be an integer 
vector. Furthermore, since the a field is not a dual gauge 
field and contains no topological excitation (otherwise the 
local constraint Eq. [Slwill be violated), the a-component 
of i for a physical topological excitation must be zero. 

As in the original quantum Hall case, The coefficients 
that appear in K and q can be understood by considering 
the EOMs resulting from it. Upon variations with respect 
to Cj, 6j, and a^ , we get: 
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(18) 

(19) 

(20) 
(21) 



The first three equations are in agreement with the pic- 
ture that spinous form integer quantum Hall states while 
holons form Laughlin v = 1/2 states under the presence 
of a flux, and that spinous carry no EM charge while 
holons carry EM charge e. Moreover, the fourth equation 
can be seen as a restatement of the occupation constraint 
Eq.[Sl 

For brevity, we shall introduce two abbreviations 
henceforth. First, we shall omit spacetime indices that 
are internally contracted. Hence we shall write eadb in- 
stead of e'^^^Gf^d^bx and [edaY instead of e'^'^^d^ax. In 
a similar spirit, we shall write dada instead of (5^0j^ — 
dua^j){dfia^ — dya^j) for the Maxwell term. Second, we 
shall write vectors and matrices in block form whenever 
appropriate, which we abbreviate by using I„ to denote 
an n-by-n identity matrix, Om,n to denote an m-hy-n 



zero matrix, and Em.„ to denote an m-hy-n matrix with 
all entries equal to 1 (such that c¥^m,n denotes an m-hy- 
n matrix with all entries equal to c). In this notation, 
the q-vector becomes q = (0;Oi^4,Oi^2,IEi,4)^ and the 
/"f-matrix in Eq. 1171 becomes: 
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( -Ei,4 
-E4.1 -I4 
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(22) 



III. 



SUPERCONDUCTING MODE AND 
PHYSICAL VORTICES 



Usually, the formation of LLs will imply that all excita- 
tions are gapped. However, this is true only if the gauge 
field is external (i.e., fixed). Since the a field is internal, 
smooth density fluctuations can occur while keeping the 
local constraint Eq. [^ and the LL structure intact. Intu- 
itively, if the a field varies across space at a sufSeiently 
long wavelength, then the spinons and holons in each lo- 
cal spatial region can still be described by the LL picture, 
but the LLs will have a larger (smaller) spacing in regions 
where the a field is stronger (weaker). Since the LL struc- 
ture is intact and the wavelength of this variation can be 
made arbitrarily long, the energy cost of such "breath- 
ing mode" can be made arbitrarily small. This breathing 
mode is thus a gapless charge-density mode of the system. 
See Fig. [3] for illustration. Note that all species of holons 
and spinons co-fluctuate with the a field in this density 
mode. A similar binding mechanism in the context of 
cuprates is proposed in Ref. UM 

The other excitations of the system can be grouped 
into two general types. The first type consists of smooth 
density fluctuations in which the fluctuations of holons, 
spinons, and a field are mismatched. The second type 
consists of quasiparticle excitations that involve holons 
or spinons excited from one LL to another. Both types 
of excitations are gapped. Since the breathing mode is 



a field 




holon 
sector 



FIG. 3: The physical picture of the breathing mode. The 
filled LL states are indicated by thick (red) horizontal lines 
while the unfilled LL states are indicated by the thin (black) 
horizontal lines. The original band structure for spinon and 
holon when no additional a. flux is also indicated in the back- 
ground (gray). 



the only gapless mode, it is non-dissipative, and hence 
the system is a superfluid when the coupling to EM fields 
are absent. Moreover, since the breathing mode includes 
the fluctuations of holons, it is charged under the EM 
field. Hence, the system will be a superconductor when 
the coupling to EM field are included^ Note that this 
superconductor breaks the time-reversal symmetry, since 
the sign of the additional amount of a flux is flipped un- 
der time reversal. Furthermore, since all four species of 
holons are binded together in the breathing mode, each 
carrying charge -he, a minimal vortex in this supercon- 
ductor is expected to carry a flux of hc/4e. We shall 
now show these claims more vigorously from the Chern- 
Simons theory Lagrangian we derived in Eq. [111 

It is easy to check that the iiT-matrix K in Eq. [IT] 
contains exactly one zero eigenvalue, with eigenvector 
Po = (2; — 2E1.4, 2Ei_2;Ei,4)"'"- Let Xi be the eigenval- 
ues of K, with Pi the corresponding eigenvectors, let 
P = [po, pi, . . . , pio] be the orthogonal matrix form by 
the eigenvectors of K, and let c' = (cg, . . . , c'^g)^ ^ ^^*^- 
Then, Eq. [TBlcan be rewritten in terms of c' as: 



+ gdc'^dc'^ + ... 

g 

= :^(q • Po)ec'„dA + {£ ■ Po)co^jy + gdc'^dc'^ + ... 

+ (terms without Cq) . 

(23) 
The Maxwell term gdc'^dc'f^ for Cq in Eq. [521 originates 
from the terms in ". . ." of Eq. I16[ which is ordinarily 
suppressed by the Chern-Simons terms. However, since 
the Chern-Simons term ecdc vanishes for Cg , the Maxwell 
terms term becomes the dominant term for Cq at low- 
energy and in the absence of external EM fields. Note 
that although the a field itself does not have a Maxwell 
term (since it arises from an occupation constraint), the 
zero-mode Cg does have a Maxwell term originated from 
the matter-field components. 

Since the Maxwell term has a gapless spectrum, we 
see that the zero-mode Cg indeed corresponds to a gap- 
less excitation. Moreover, since all other gauge-field com- 



ponents have non-zero Chern-Simons terms, excitations 
in these gauge-field components are gapped (these exci- 
tations corresponds to the "mismatched" density fluctu- 
ation mentioned earlier), verifying the earlier assertion 
that there is only one gapless density mode. Moreover, 
since q • po 7^ 0, we see that the zero-mode is indeed 
charged under the external EM field. Hence, as argued 
above, the doped system is a superconductori^i 

The eigenvector pi can be interpreted as the ratio of 
density fluctuations between the different field compo- 
nents in the mode c[ . Thus the zero-mode indeed involves 
the fluctuations of all species of spinous and holons, tied 
together by the internal a field. 

Since the system is a superconductor, when a suffi- 
ciently large external B field is applied, physical vortices, 
with the amount of flux through each vortex quantized, 
are expected to form. In the Chern-Simons formulation, 
these physical vortices manifest in the topological term 
(i.e.; the {i ■ c)^ijy term) in Eq. [TH Taking an isolated 
topological excitation with (jy,jy,jy) = ((5(a;— xq), 0,0), 
considering the EOM associated with Cg as resulted from 
Eq. [231 and remembering that {edAf = e^^'''^f,A^ = B 
is the physical magnetic field, we obtain (in units which 
h = c=l): 
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5{x ~ Xg) -f 



(24) 



This is the Meissner effect, which again confirms that 
the system is a superconductor. Moreover, it is easy to 
check that non-zero |(^-po)/(q-Po)| has a minimum of 1/4 
(attained by, e.g., an .^-vector having a single "-I-1" in one 
of its bj components and "0" in all its other components). 
From this we conclude that the magnetic flux through a 
minimal vortex is hc/Ae, justifying the intuitive claim 
given above. 



IV. QUASIPARTICLES — STATISTICS 

It is important to note that not all topological exci- 
tations are EM-charged. The structure of these EM- 
ncutral topological excitations highlights the differences 
between this system and a conventional superconduc- 
tor, and hence qualify the adjective "exotic." We shall 
call these EM-neutral topological excitations "quasiparti- 
cles," to distinguish them from the EM-charged "physical 
vortices" considered in the previous section. 

From Eq. [Ml a topological excitation carries a non- 
zero magnetic flux if and only if £ • po 7^ 0. In other 
words, a topological excitation is EM-neutral if and only 
if it docs not couple to the zero-mode. Note that quan- 
tity (. ■ pg can be regarded as the zero-mode "charge" 
carried by the topological excitation. A topological ex- 
citation with .^ • Pg 7^ couples to the zero-mode and 
carries its "charge," which induces an 1/r "electric" field 
of the zero-mode and gives rise to a diverging energy gap 
A ~ Ini, where L is the system size. In comparison, a 
topological excitation that satisfies ^ • po = is decoupled 



from the zero-mode and hence has a finite energy gap and 
short ranged interactions. These EM-neutral topological 
excitations are thus analogous to the (possibly fraction- 
alized) quasiparticles in quantum Hall systems, and it is 
sensible to consider the (mutual) statistics between them. 
Recall that the set of £- vectors (which may have non- 
zero a-component) form an eleven dimensional vector 
space. The set of ^-vectors satisfying ^-po — forms a ten 
dimensional subspace of this eleven dimensional space. 
The iiT-matrix restricted to this subspace, K,., is invert- 
ible. Hence we can integrate out the gauge fields associ- 
ated with this subspace (i.e., the gauge fields c[, . . . , c'^q 



in Eq. 



This will convert the terms we omitted in 



Eq. [23] under the texts "terms without Cq" into a Hopf 
term. Explicitly, upon integrating out c'l, . . . , c'^q the La- 
grangian takes the form: 

>C" = :f (q • Po)ec'odA + {i ■ po)c^^.?(; + gdc'odc'o + ... 

Ztt 

+ .(FrA--(^)? + ... 

= (terms with <:[,) + ^(F)^^^' (^^) J' + ' ' ' ' 

(25) 
(c.f. Eq.[TlD, where j'^ = j!^e + {e/2Tr)iedAyq. 

As in the quantum Hall case, from Eq. [21] the statis- 
tical phase 9 when one quasiparticle described by jyi 
winds around another described by j'y^l' can be read 
off as = 2tt £'^ K~^£' . For identical quasiparticles, 0/2 
gives the statistical phase when two such quasiparticles 
are exchanged. 

For explicit computation a basis for £- vectors for this 
ten-dimensional subspace must be specified. Naively one 
may simply choose this basis to be the set of eigenvectors 
of K having non-zero eigenvalues. This choice turns out 
to be inconvenient as some of the eigenvectors of K are 
non-integer while the quantization condition requires all 
£ to be integer vectors. Hence, instead we shall use the 
following basis: 



h = {0 

f2 = (0 
f3 = (0 
?4=(0 
^5^(0 
^ = {0 
f7=(0 
f8 = (0 
?9 = (0 



-l,l,0,0,Oi,2;Oi,4)^ 
-1,0,1,0,01,2; 01,4)^ 
-l,0,0,l,Oi,2;Oi,4)^ 
Oi.4,Oi,2;0,0,l,-l)^ 
O1.4, 01,2; 0,1, 0,-1)^ 
Oi.4,Oi,2;l, 0,0,-1)^ 
0,l,0,0,Oi,2;0,l,l,0)^ 
1,0, 0,0, 1,0; 01,4)^, 



(26) 



=-1,4 



IE1.2 



hAf, 



^io = (-l;Oi,4,0,l; 



"lA) 



It can be shown that all integer ^-vectors satisfying 
^ • Po = can be written as integer combinations of 
the above basis vectors. It should be remarked that £1 



through £q are indeed eigenvectors of K, with £1 through 
£3 having eigenvalue — 1 and £4 through £q having eigen- 
value 2. However, ^7 through ^10 are not eigenvectors of 
K. 

In this basis, Kl^^ takes the form: 



k: 



/ -2 -1 -1 


-1 


1 \ 


-1 -2 -1 O3.3 





1 ©3,2 


-1 -1 -2 
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1 1/2 1/2 
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©3,3 1/2 1 1/2 
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©3,2 


1/2 1/2 1 
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©2,3 ©2,3 






1 
11/ 



(27) 

Note that ^10 contains a non-zero a-component and is 
thus unphysical. Moreover, from our interpretation of as 
and ag as arising from the physics of band bottoms, we 
expect a topological excitation in these two components 
to be much more energetically costly than those of the 
other matter fields. Hence we can also neglect £s and £q. 
Thus only the top- left block of K~^ is relevant for the 
statistics of low-lying physical quasiparticle excitations. 
Henceforth we shall restrict the meaning "quasiparticle" 
to those whose ^-vector is an integer combination of £1 
through £t. 

From K~^ it can be seen that the system contains 
quasiparticles with non-trivial mutual statistics. In par- 
ticular, there there are fermions having semionic mutual 
statistics (i.e., a phase factor of tt when one quasiparticle 
winds around another), manifesting in, e.g., quasiparti- 
cles described by £4 and £c, . 

The self-statistics and mutual statistics of different 
quasiparticles can be understood intuitively. Recall that 
our system is constructed by coupling integer and frac- 
tional quantum Hall states via a common constraint 
gauge field a. If we assume that the different quantum 
Hall states are independent of each other, i.e., a "charge" 
in one matter-field component has trivial bosonic statis- 
tics with a "charge" in a different matter-field compo- 
nent, then the statistics of these quasiparticles can be 
read off by considering their underlying constituents. For 
example, since ^4 and £5 overlaps in one v = 1/2 compo- 
nent, their mutual statistics is semionic. Similarly, since 
£4 overlaps with itself in two v ^ \/2 components, its 
self-statistics is fermionicr^S From this intuitive picture, 
it is evident that a "+1" in a spinon component in the 
^-vector should be identified with a spinon excitation on 
top of the integer quantum Hall state that formed near 
the corresponding Dirac node, while a "-|-1" in a holon 
component in the -^-vector should be identified with half- 
holon excitation on top of the v = 1/2 quantum Hall 
state that formed near the corresponding band bottom. 
Similarly, a "—1" in a spinon (holon) component in the 
^-vector should be identified as an anti-spinon (anti-half- 



lialf-holon half-anti-holon 




^ 




(a) 



anti-spinon 
(b) 



holon anti holon 

(c) (d) 



FIG. 4: Physical interpretation of ^-vector: (a) a single "+1" 
in a spinon component identified as spinon; (b) a single "—1" 
in a spinon component identified as anti-spinon; (c) a single 
"+f" ("+2") in a holon component identified as half-holon 
(holon); and (d) a single "—1" ("—2") in a holon component 
identified as anti-half- holon (anti-holon). The thick (red) hor- 
izontal lines indicate filled LLs that forms the ground state 
of the system, while the thin (black) horizontal lines indicate 
unfilled LLs. 
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FIG. 5: Visualization of the (a) "elementary" SQP; (b) "el- 
ementary" HQP; (c) "minimal" MQP of the first type; and 
(d) "minimal" MQP of the second type. 



holon). Sec Fig. m for illustration. 

To discuss these quasiparticlcs further, it is useful to 
divide them into three classes. The first class consists 
of quasiparticles with spinon components only and will 
be referred to as "spinon quasiparticles" (SQP). The sec- 
ond class consists of quasiparticles with holon compo- 
nents only and will be referred to as "holon quasipar- 
ticles" (HQP). The remaining class consists of quasi- 
particles that have both spinon and holon components, 
and will be referred to as "mixed quasiparticles" (MQP). 
The first two classes can be constructed by compound- 
ing "elementary" quasiparticles of the same type. For 
SQP, the "elementary" quasiparticles are described by 
^-vectors having exactly one "-I-1" component and one 
"—1" component in the spinon sector (e.g., the £i,i2, 
and £3 in Eq. \7^ . For HQP, the "elementary" quasipar- 
ticlcs are described by £- vectors having exactly one "-I-1" 
component and one "—1" component in the holon sec- 
tor (e.g., the ^4,^5, and 4 in Eq. ^. As for the MQP, 
one can start with "minimal" quasiparticles with exactly 
one "-I-1" component in the spinon sector and one "+2" 
components in the holon sector, and build all MQP by 
compounding at least one such "minimal" quasiparticles 
together with zero or more "elementary" SQP and HQP. 
Alternatively, one may start with a second type of "mini- 
mal" quasiparticle in the MQP sector, which has exactly 
one "-I-1" component in the spinon sector and two "-I-1" 
components in the holon sector, and build all MQP by 
compounding at least one such "minimal" quasiparticles 
together with zero or more "elementary" SQP and HQP 
(note that the second-type of "minimal" MQP is simply 
a "minimal" MQP of the first type compounded with 
an "elementary" HQP. The introduction of two different 
types of "minimal" MQP will be clear in the following). 

These "elementary" and "minimal" quasiparticle ex- 
citations can be visualized in the following way: The 
"elementary" SQP can be visualized as a particle-hole 
excitation in the spinon quantum Hall levels, in which 
a spinon is removed from one Dirac node and added in 



TABLE I: Self- and mutual- statistics of the "elementary" or 
"minimal" quasiparticles in the doped kagome system. The 
adjective "elementary" or "minimal" are omitted but assumed 
in the table entries. The subscript I and II indicates the type 
of "minimal" MQP considered (see the main text for their 
definitions). When an entry contain multiple cases, both cases 
are possible but are realized by different quasiparticles in the 
respective sectors. 



Type 


Self-Statistics" 


Mutual Statistical Phase' 
SPH HPH MQPi MQPii 


SQP 


b 


2tt 


2-K 


2-K 


2-K 


HQP 


f 


2-K 


-K or 2tt 


2-K 


-K or 27r 


MQPi 


f 


27r 


2n 


2-K 


27r 


MQPii 


b 


2-JT 


■K or 2-K 


2-K 


-K or 27r 



''b=bosonic, f=fermionic, s=semionic 

'Phase angle accumulated when one quasiparticle winds around 
another, modulo 2-k. 



another. The elementary HQP can be visualized as a 
particle-hole excitation in the holon quantum Hall lev- 
els, in which a half holon is transferred from one band 
bottom to another. The minimal SQP can be visualized 
as adding both spinon and (half) holons into the original 
system. See Fig. [5] for illustrations. 

With this classification, the information on the self- 
and mutual- statistics of the quasiparticlcs contained in 
K~^ can be summarized more transparently in terms of 
the self- and mutual- statistics of the "elementary" SQP, 
"elementary" HQP, and "minimal" MQP. The result is 
presented in Tabled 



V. QUASIPARTICLES— QUANTUM NUMBERS 

Since the quasiparticlcs have finite energy gaps and 
short-ranged interactions, they may carry well-defined 
quantum numbers. In particular, it is sensible to con- 



sider the k quantum numbers for these quasiparticles, 
since they arise from LLs that form near Dirac points 
or band bottoms with well-defined crystal momentum k. 
Similarly, it is sensible to consider the Sz quantum num- 
bers for quasiparticles with spinon components. We shall 
see that this program can be carried out for "elementary" 
spinon quasiparticles and for the "minimal" mixed quasi- 
particles of first type, but not easily for the "elementary" 
holon quasiparticles and the "minimal" mixed quasipar- 
ticles of the second type. 



Recall that we constructed a tight-binding model with 
doubled unit cell for the DSL ansatz. The unit cell is nec- 
essarily doubled because the DSL ansatz enclose a flux 
of TT within the original unit cell spanned by ri/2 — x 
and Y2 = (l/2)x + (\/3/2)y (c.f. Fig. |l(a)[), and hence 



the operators that corresponds to translation by x, Tx, 
and the operator that corresponds to translation by r2, 
Tr2, do not commute in general (i.e., \Tx, T'r2] 7^ 0), even 
though both commute with the mean-field tight-binding 
Hamiltonian. Consequently, single-spinon and singlc- 
holon states in the DSL ansatz generally form multi- 
dimensional irreducible representations under the joint 
action of T^ and Tr2 (i-s., T^ and Tr2 manifest as multi- 
dimensional matrices that cannot be simultaneously di- 
agonalized when acting on these states), and cannot be 
labeled simply by a pair of numbers (ci, C2) as in the or- 
dinary casci^ Furthermore, the matrices for T^ and Tr2 
will in general be a-gauge-dependent. However, when an 
even number of spinon and holon excitations are consid- 
ered as a whole, the total phase accumulated when the 
particles circle around the original unit cell becomes a 
multiple of 27r, and thus [Tx,Tr2] = m such subspace. 
Hence it is possible to reconstruct the crystal momentum 
in the original Brillouin zone if our attention is restricted 
to such states. The tool for reconstructing the crystal 
momentum in the original Brillouin zone is known as 
the projective symmetry group (PSG)J^ Physically, the 
gauge dependence of single-spinon and single-holon states 
indicate that they cannot be created alone. 



It can be checked that all SQP are composed of an even 
number of spinous and anti-spinons. The above discus- 
sion then implies that they carry well-defined k quantum 
numbers in the original Brillouin zone. To derive the 
transformational properties under T^ and Tr2, we com- 
pute the transformation properties of the original spinon 
matter fields. The procedures for doing so have been de- 
scribed in details in Ref. 0, here we shall just state the 
results. 



Let 771 , . . . , 774 denote the topological excitations near 
the four (two k-vectors and two spins) Dirac nodes as 
indicated in Fig. 



6(a) Then, assuming that they have 
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FIG. 6: (a) The labels for the four spinon topological excita- 
tions, (b) Physical interpretation of the "missing states" in 
the fixed S~ quantization and doubled unit cell Chern-Simons 
formulation, (c) Spectrum of "elementary" SQP, with k and 
Sz quantum number indicated, before restoring full symme- 
try, (d) Spectrum of "elementary" SQP after restoring the 
SU(2) symmetry by adding extra quasiparticles. The dotted 
arrows indicate equivalent k-point upon translation by the 
original reciprocal lattice vectors (spanned by 2ki and k2 in 
Fig. 1(b) I. For dimension of the Brillouin zone, c.f. Fig. 1(b) 



spinon fields at the same Dirac nodes. 



TAm] 


= eWi2^,, 


Tr2[ril] 


^e"lW. 


TAm] 


= e"-/i2^i. 


Tr2[m] 


- e— /^/2 


T.[m] 


= eWi2,,,, 


Tr2[m] 


==e'-/^%, 


TxVia] 


= e"Wi2^3^ 


Tr2[m] 


= e--/2r;4 



(28) 



Furthermore, we assume that T.^ and Tr2 satisfy the 
generic conjugation and composition laws: 



T{r\ = {Tm* , T{^ ■ ^'] = T{^A ■ TW\ 



(29) 



where t/), "i/)' denotes generic quasiparticle states, "0* de- 
notes an anti-particle of '0, and V' ' V'' denotes a bound 
state composed of ^ and V''- 

A general basis for "elementary" SQP is spanned by 
r\ir\* with i ^ j ■ There are twelve distinct "elementary" 
SQP, which form six reducible representations under T^ 
and Tj.2- Upon diagonalization, the resulting "elemen- 
tary" SQP in the new basis each carry distinct Sz and k 
(in the original Brillouin zone) quantum numbers. These 



are summarized in Fig. 6(c) 



the same transformational properties as the underlying 



Notice that Fig. 6(c) is somewhat unsettling. First, 
even though we have not performed a PSG study on 
rotation operators, intuition on rotation symmetry sug- 
gests that there should be four states (with Sz = 



10 



1, — 1, 0,and 0) located at k = (tt, — 7r/-\/3). Second, al- 
though our Chern-Simons theory is formulated with a 
fixed quantization axis for spin, the SU{2) spin-rotation 
symmetry should remain unbroken. Therefore, the Sz 
eigenvalues should organize into SU{2) representations 
for each k value. While this is true for k = (tt, 7r/\/3) 
and k — (0, 27r/V3); where the "elementary" SQP form 
10 representations, the same does not hold for k = 
(7r,-7r/\/3) andk= (0,0). 

The two issues mentioned above indicate that some 
topological excitations arc lost in our formulation. In 
other words, there are topological excitations that have 
trivial Sz quantum numbers but non-trivial S quantum 
numbers. Similarly, there are topological excitations that 
have trivial k quantum numbers in the reduced Brillouin 
zone but non-trivial k quantum numbers in the original 
Brillouin zone. Physically, the original of these missing 
excitations can be understood as follows: in the hydro- 
dynamic approach, an £- vector with a single "-I-1" in a 
spinon component represent a spinon at a Dirac node, 
while £- vector with a single "—1" in a spinon component 
represent an anti-spinon at a Dirac node. The previously 
defined set of ^-vectors that characterized the "elemen- 
tary" SQP fail to captured an excitonic state in which a 
spinon is excited from a filled LL to an empty LL, thus 
leaving an anti-spinon behind (see Fig. |6(b)| for an illus- 
tration), which precisely carry trivial Sz quantum num- 
bers and transoform trivially under T^i and Tr2- Note 
that there are four possible excitonic states of this form, 
hence we expect four states to be added. In our Chern- 
Simons formulation, these excitations may be disguised 
as combinations of density operators (~ 9c). 

From Fig. 6(c) and the forgoing discussions, it is evi- 
dent that extra states should be added at k = (tt, 7r/\/3) 
and k = (0,0), so that the states at k = (0,0) and 
k = (tt, tt/^/3) each form a 1©0 representation of SU{2). 
The final result after making this reparation is shown in 
Fig. [DJc). Formally, the same result can be reached if 
we allow objects of the form rjiTj* to be counted as ele- 
mentary SQP, then apply Eq. [^5] and the procedure of 
diagonalization as before in this extended basis. 

Observe that the "elementary" spinon SQP (and hence 
the entire SQP sector) all carry integer spins. However, 
we also know that a conventional superconductor con- 
tains spin-1/2 fermionic excitations (i.e., the Bogoliubov 
quasiparticles). From our assignment of Sz quantum 
number and from the table of quasiparticle statistics Ta- 
ble [H it is evident that the "minimal" MQP of the first 
type play the role the these Bogoliubov quasiparticles in 
the doped kagomc system. In contrast, minimal MQP 
of the second type are spin-1/2 quasiparticles that carry 
bosonic statistics and hence is another distinctive signa- 
tures of this exotic superconductor. 

Since a "minimal" MQP of the first type can be treated 
as a bound state of a spinon and a holon (c.f. Fig. 5(c) I, 
the k quantum number in the original Brillouin zone are 
again well-defined for them. To construct their quantum 
numbers, we need to know how holons transform under 




FIG. 7: (a) The labels for the four holon excitations, (b) 
Spectrum of "elementary" MQP, with k quantum number in- 
dicated. Each point in k space forms a S = 1/2 representation 
in spin. The dotted arrows indicate equivalent k-point upon 
translation by the original reciprocal lattice vectors. 



Tx and Tr2- Let (/5i, . . . , (^4 denotes the half-holon exci- 
tations near the four holon band bottom as indicated in 
Fig. 6(a) such that (pi, . . . ,if>\ denotes the correspond- 



tM] - 


= ^L 


T.2[^?] = e-/6^?, 


tM] 


= V3. 


T.2bi]-e— /VL 


tM] -- 


= e-'-'W , 


T^M]-e''^/'^l, 


tM] -- 


-- e'"/V? , 


r.2[^I] = e-wV^ 



ing holon excitations (c.f. Fig. 4(c)). Following the same 
procedure that produces Eq. 1281 we obtain the transfor- 
mation laws: 



(30) 



A general basis for "minimal" MQP of the first type 
is spanned by rjiip^,. There are sixteen distinct first- type 
"minimal" MQP, which form eight reducible representa- 
tions under T^ and Tr2- Upon diagonalization, the re- 
sulting first-type "minimal" MQP in the new basis each 
carry distinct Sz and k (in the original Brillouin zone) 
quantum numbers, and the full SU(2) representation in 
spin can be recovered trivially by combining spin-up and 
spin-down states. The final results are summarized in 
Fig. [7(b)] 

Having considered the SQP sector and the "minimal" 
MQP of the first type, one may attempt to carry out sim- 
ilar analysis for the HQP sector and for the "minimal" 
MQP of the second type. However, in doing so, issues 
arise from the fractionalization of holons into half- holons. 
Recall that in deriving the transformational rules of the 
quasiparticles, we identify the components of i as be- 
ing spinon and holon excitations, and assume that these 
excitations carry the same quantum numbers as the un- 
derlying spinous and holons that form the LLs in the first 
place. However, the HQP sector and the "minimal" MQP 
of the second type are bound states that involve half- 
holons, whose quantum numbers cannot be directly in- 
ferred from the underlying spinons and holons. More con- 
cretely, we need to know the transformation laws T[ifiLp*] 
for half-holon-anti-half- holon pairs ^Pi^p* in order to con- 
struct their quantum numbers, but we only have informa- 



11 



tion about transformation laws ^^[(/j^] of holon excitation 



VI. AN ALTERNATIVE DERIVATION BY 
ELIMINATING THE AUXILIARY FIELD 



It is far from clear how TYpiLp*] can be related to T[93^]. 
The answer for such question may even be non-unique. 
We have already seen an analogous situation in the for- 
going discussion: while the spinon-anti-spinon pairs rjiT]* 
have well-defined gauge-invariant k quantum numbers in 
the original Brillouin zone, the single spinous 77,; form 
gauge-dependent two-dimensional representations under 



The possible ambiguity in the transformation law 
T[LpiLp*] of half-holon-anti-half-holon pairs V'iV^ signifies 
that it may not be possible to produce these quasipar- 
ticles alone. Although a half-holon-anti-half-holon pair 
can be thought of as resulted from removing a half-holon 
from one band bottom and adding one in another, it 
is not clear that the process can be done in via single 
half-holon tunneling. This is analogous to the case when 
two fractional quantum Hall system are separated by a 
constriction, where it is only possible to tunnel physical 
electrons ji^ 



Combining the results from Sect. IIVI and |Vl we see 
that there are two very different class of quasiparticle 
excitations in the doped kagome system — which can be 
termed as "conventional" and "exotic," respectively. The 
"conventional" class consists of quasiparticles that can 
be created alone, which carry well-defined crystal mo- 
mentum k in the original Brillouin zone and possess con- 
ventional (fermionic or bosonic) statistics. These include 
the spinon particle-holes, the holon (but not half-holon) 
particle-holes, the "minimal" mixed quasiparticles of the 
first type (a.k.a. the "Bogoliubov quasiparticles"), and 
their composites. In contrast, the "exotic" class con- 
sists of quasiparticle that cannot be created alone, whose 
crystal momentum may not be well-defined, and whose 
statistics may be fractional. These include the half-holon 
particle-holes and the "minimal" mixed quasiparticles of 
the second type (which are "Bogoliubov quasiparticles" 
dressed with a half-holon particle-hole). In terms of the 
underlying electronic system, the former class are excita- 
tions that are local in terms of the underlying electron 
operators c and c^, while the latter class are excitations 
that are non-local in terms of c and c^. 



It should be warned that questions regarding the ener- 
getics (and hence stability) of the quasiparticles have not 
been touched in Sect. IIVI and FVl In particular, it is not 
clear whether the bosonic or the fermionic spin-1/2 exci- 
tation has a lower energy. Though this information is in 
principle contained in the Maxwell term Eg. 1101 to obtain 
it requires a detailed consideration of the short-distance 
physics in the t-J model, and is beyond the scope of this 
paper. 



It is a curious result that in Eq. 1271 once the unphys- 
ical ^10 is removed from the spectrum, the quasiparticle 
represented by ig becomes purely bosonic (i.e., having 
trivial bosonic mutual statistics with all other quasipar- 
ticles and trivial bosonic self-statistics). This suggests 
that ^9 corresponds to some local density excitation of 
the system and thus should not be regarded as topologi- 
cal. Moreover, the procedure of first treating ^10 as part 
of the spectrum in computing K^^ and then removing 
this degree of freedom at the very end of the calculation 
seems somewhat dubious. Recall that the gauge field 
a^ is introduced to enforce the occupation constraint 
Eq. [51 This gauge field is thus an auxiliary field that 
is void of self-dynamics (i.e., the term eada vanishes) 
and topologically trivial (i.e., the a-componcnt of ^ must 
be zero). Therefore, one may attempt to re- derive the 
previous results by eliminating this a field right at the 
beginning by enforcing the constraint directly. This can 
indeed be done, as we shall show in the following. 

Recall that the EOM with respect to a^ leads to the 
constraint equation Eq. [5T] in the Chern-Simons formu- 
lation. From this, one may argue that the effect of intro- 
ducing the a field can alternatively be produced by setting 
^j flj +X]j &/ == directly. To do so, we perform a two- 
step transformation on the Lagrangian Eg. 1161 First, we 
set: 



a^^ = XI "/ ^n^ "'/ = «/ f°i' ^ ^ 6 ; (31) 
fe'i^ = X «J and h'/ = 6^ for J 7^ 1 . (32) 



Then the constraint becomes Og^ 
enforce directly by setting: 



= 6'/ 



h'l^ = 0, which we 



(33) 



thus eliminating one variable. 

Note that since the a field appears in Eq. [1^] only 
through the term e{J2i^i + J2j^j)^'^j i^ 8°^ dropped 
out of the transformed Chern-Simons Lagrangian. Let- 
ting c — {p; a'l, . . . , flg; 621 • • • J ^4)7 which is a column vec- 
tor of only nine (as opposed to eleven) gauge fields. Eg. [TCI 
becomes: 



C = -Lec^A-ae - ^e(q • c)dA + (i ■ c)^jf^ + ..., (34) 

47r ZTT 



where q — (1; O1.5; ©1^3)^ is the transformed charge vec- 



12 



tor, and K is the transformed JC-matrix: 



K = 



( ^ 


11111 


-2 


-2 -2\ 




1111 








10 111 








110 11 




©5,3 




1110 1 








11112 






-2 




4 


2 2 


-2 


©3,5 


2 


4 2 


\-2 




2 


2 4 / 



(35) 



As for topological excitations, from the transformation 
between c and c, it can be seen that the correspondence 
between i and i reads: 



(0;nai, . . .na6;nbi, . . .nu) 



= inbi-nae;nai-na6, ■ ■ • jWaS-'^ae; 



(36) 



Hence, (. is an integer vector if and only if (. is also an 
integer vector. Moreover, from Eq.[3niit can be seen that 
£9 is mapped to ^ = 0^^ which is consistent with our 
previous argument that the quasiparticle corresponding 
to ^9 is purely bosonic and hence should be considered 
as non-topological. 

Although K and K look rather different superficiallyf^^ 
all the major conclusions from Sect. IIIlHVl can be repro- 
duced with K. In particular, we shall check that the ex- 
istence of a single gapless mode, the hc/Ae flux through a 
minimal vortex, and the semionic quasiparticle statistics 
can all be obtained from K . 

It is easy to check that K has exactly one zero eigen- 
value, with po ~ (4; — 2Ei 4, 2; £13)-^ its eigenvector. Us- 



ing the transformation equations Eq. [3TH551 we see that 
this eigenvector corresponds precisely to the eigenvector 
Po we found in Sect. IIIII Thus, again we conclude that 
the system contains a gapless mode associated with su- 
perconductivity, and that this gapless mode can be inter- 
preted as fluctuations of all spinous and holons species 
whose ratio is matched (through their common coupling 
to the gauge field a^^). 

Moreover, the amount of magnetic flux that passes 
through a physical vortex is still described by Eq. [24] 
upon the obvious modifications. Since q • po = 4, we 
recover the conclusion that a minimal physical vortex 
carries a flux of hc/Ae. Furthermore, it can be checked 
that ^ • Po = ^ • Po for £, (. satisfying the correspondence 
Eq. [351 Hence the flux carried by a vortex calculated 
from K agrees with the value calculated from K. 

As before, the quasiparticle excitations (which are EM- 
neutral, short-ranged interacting, and have finite energy 
gaps) are characterized by the condition that ^ • po = 0, 
which defines an eight-dimensional subspace of the nine- 
dimensional space in this case. The A'-matrix restricted 



to this subspace, AT^, is invertible. We may choose a basis 
for this subspace that corresponds to the basis choice 
Eq. [2n|in the original representation. Explicitly, 



K-^ = 



li = (0;0, -1,1,0,0; 


Oi,3)^ 






£2 = (0; 0,-1, 0,1,0; 


01,3)^ 






^3 = (0;0, -1,0,0,1; 


Oi,3)^ 






^4 = (0; 0,01,4; 0,1, - 


-1)^ 






4 = (0; 0,01,4; 1,0, - 


-1)^ 






^6 = (1; 0,01,4; -1,- 


1,-2) 


T 




£7 = (0;0,0,1,0,0;1, 


1,0)^ 






4 = (0;l,l,0,0,0;O 


1,3)^ 






:;an be checked that: 








-2 -1 -1 




-1 


l\ 


-1 -2 -1 03,3 







1 


-1 -1 -2 







1 


1 1/2 


1/2 


1/2 





03,3 1/2 1 


1/2 


1/2 





1/2 1/2 


1 








-10 1/2 1/2 











1110 








0/ 



(37) 



(38) 



in agreement with the results in Sect. ITVl 



VII. CONCLUSIONS 

In this paper we have considered the theory of a doped 
spin-1/2 kagome lattice described by the t-J model. We 
start with the slave-boson theory and the assumption 
that the undoped system is described by the U{1) Dirac 
spin liquid, from which we argued that the doped system 
is analogous to a coupled quantum Hall system, with 
the role of the external magnetic field in the usual case 
taken up by an emergent gauge field a. The analogy 
with quantum Hall systems compels us to introduce the 
Chern-Simons theory as an effective description of the 
low-energy physics of the system. This allows us to de- 
scribe the superconductivity, the physical vortices, and 
the electromagnetically neutral quasiparticles in a unified 
mathematical framework. We show that there are two al- 
ternative Chern-Simons theories that produce identical 
results — one with the auxiliary field a kept until the end, 
and the other with the auxiliary field and a redundant 
dual matter field eliminated at the beginning. 

In our scenario, the coupled quantum Hall system con- 
sists of four species of spinous and four species of holons 
at low energy. We show that such system exhibit super- 
conductivity and that the flux carried by a minimal vor- 
tex is hc/Ae. The system also contains fermionic quasi- 
particles with semionic mutual statistics, and bosonic 
spin-1/2 quasiparticle. As for the quantum numbers car- 
ried by the quasiparticles, we analyzed the spinon sector 



13 



in details and found that it is possible to recover the full 
SU{2) and (un-enlarged) lattice symmetry of the "ele- 
mentary" quasiparticles in this sector, upon the inclu- 
sion of quasiparticles that are not easily represented in 
the original fixed-spin-quantization-axis, enlarged- unit- 
cell description. The same classification of quantum num- 
bers are also carried out for the spin- 1/2 fcrmionic quasi- 
particles, which are the analog of Bogoliubov quasiparti- 
cles in our exotic superconductor. 

In this paper we have argued that the doped spin- 
1/2 kagome system may exhibit exotic superconductivity 
that is higher unconventional. However, it should be re- 
marked we have presented only one possible scenario for 
the doped kagome system. For example, it is possible 
that the ground state of the undoped system is a valence 
bond solidS, and hence invalidate our analysis. Further- 
more, experimentally realizing the idealized system con- 



sidered considered in this paper may involve considerable 
difhculties. For instance, in the case of Herbertsmithite, 
it is known that the substitution between Cu and Zn 
atoms can be as big as 5% Ji It is our hope that this pa- 
per will generate further interests in the doped spin-1/2 
kagome system, as well as other systems that may ex- 
hibit anlogous exotic superconducting machanisms, both 
experimentally and theoretically. 
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Anderson-Higgs mechanism in which this Goldstone mode 
is "eaten up" by the electromagnetic field. 

^^ The existence of the zero-mode (and hence superconduc- 
tivity) is in fact a rather general consequence of zero total 
HaU number (i.e., Eaii species " = 0)- See ReiM- 

^^ For this intuitive picture to be accurate, the sign of the 
component must also be taken into account. 
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and Tr2, respectively, and are related to the crystal mo- 
mentum k in the original Brillouin zone via exp(ik-£) — c\ 
and exp(zk ■ r2) = C2. 

^* More generally, given (. in the transformed basis, the cor- 
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■^^ It can even be checked that K contains irrational eigenval- 
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